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Abstract. We compute the number level variance £2 and the level compressibility \ from high precision 
data for the Anderson model of localization and show that they can be used in order to estimate the 
critical properties at the metal-insulator transition by means of finite-size scaling. With N, W, and L 
denoting, respectively, linear system size, disorder strength, and the average number of levels in units of 
the mean level spacing, we find that both x(N, W) and the integrated £2 obey finite-size scaling. The high 
precision data was obtained for an anisotropic three-dimensional Anderson model with disorder given by 
a box distribution of width W/2. We compute the critical exponent as v ~ 1.45 ± 0.12 and the critical 
disorder as W c ~ 8.59 ± 0.05 in agreement with previous transfer-matrix studies in the anisotropic model. 
Furthermore, we find \ ~ 0.28±0.06 at the metal-insulator transition in very close agreement with previous 
results. 
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1 Introduction 



The Anderson metal-insulator transition (MIT) in disor- 
dered systems has been vigorously studied for a long time 
and still continues to attract much attention jj] . For 
non-interacting electrons in disordered systems the scaling 
hypothesis of localization has been successfully validated 
by theoretical ||[|] and numerical [fz| — p~0|| approaches. The 
latter approaches use well-known techniques of finite-size 
scaling (FSS) Q. FSS at the Anderson MIT has a note- 
worthy history, reaching a first peak with the seminal 
papers of Pichard/Sarma iR Bl and MacKinnon/Kramer 
R| . Especially in Ref . (Rj , the groundwork for a reli- 
able, numerical FSS procedure was laid and scaling curves 
could be constructed that proved the existence of an MIT 
in 3D and the absence of such in 2D and ID. In these and 
later studies based on the same analysis technique Q] , the 
critical exponent as estimated from the divergence of 
the infinite-size localization (correlation) length £(W) as 
a function of the disorder strength W at the transition 
W = W c , i.e., £ oc |1 — W/W c \~ v , was systematically un- 
derestimated. The divergent nature at the transition could 
only be poorly captured by FSS of data obtained for small 
system sizes and large errors e in these finite-size data. 
However, as more powerful computers became available 



in the last decade, one observed a trend towards larger 
values of v « R35 for e < 1%. 

In 1994, high-precision data (e < 0.2%) of MacK- 
innon |r| for the Anderson model of localization (AM) 
showed a hitherto neglected systematic shift of the transi- 
tion point W c with increasing system size. Taking this into 
account phcnomenologically, v — 1.54 ± 0.08 was found 
JL6| . A subsequent approach by Slevin/Ohtsuki [p7|-pj[ 
incorporated these shifts as a consequence of irrelevant 
scaling variables and further allowed for corrections to 
scaling due to nonlincarities. With higher-precision data 
(e w 0.1%), they found v = 1.57 ± 0.04. Further re- 
sults for, e.g., the AM with anisotropic hopping [p0|-p2|, 
the off-diagonal AM ||||24|], and the AM in a magnetic 
field [2li,Bq], confirmed this value of v within the error 
bars [27]. Also, v is identical for the MIT as a function 
of disorder or energy E3, 24 1. We emphasize that a prop- 
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erly performed Slevin/Ohtsuki scaling procedure needs to 
assume various fit functions and that the final estimates 
are to be suitably extracted from many such functional 
forms and starting parameters pl|-]23]] ; bootstrap p7[-^9[ 
or Monte Carlo methods |fil|-f23"f then need to be employed 
for a precise estimate of error bars. 

Regarding experiments, we note that similarly precise 
data (0.1%) are much harder to obtain for our experimen- 
tal colleagues. Nevertheless, recent advances in this direc- 
tion based on careful finite-temperature analysis of the 
conductivity data show a clear trend towards increasing 
v > 1 B2j. The roles of sample inhomogeneities, mag- 
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netic effects and other possible experimental influences are 
also discussed 

The statistical properties of spectra of disordered single- 
electron systems are closely related to the localization 
properties of the corresponding wave functions p6|-|38|]. 
In the 3D AM we have the insulating, the critical and the 
metallic phases, respectively. For the insulating regime, 
the localized states even if they are close in energy have 
an exponentially small overlap and their levels are uncor- 
rected. Accordingly, in the thermodynamic limit the nor- 
malized distribution of the spacing s between neighboring 
energy levels follows the Poisson law 



Pp(s) = exp(-s). 



(1) 



In the metallic regime, the large overlap of delocalized 
states induces correlations in the spectrum leading to level 
repulsion. In this case, if the system is invariant under ro- 
tational and under time-reversal symmetry, the normal- 
ized spacing distribution closely follows the Wigner sur- 
mise of the Gaussian orthogonal ensemble (GOE) of ran- 
dom matrices 



Avd(s) = 77 s ex P 



"4 S 



(2) 



The third symmetry class at the MIT is usually called the 
critical statistics p5| , p^ -|4^|. Its normalized level spacing 
distribution for large s is proportional to 



P c (s) oc exp (— ks) 



(3) 



where k has been argued to be related to the value of the 
level compressibility Xc at the MIT jf7| . 

Various measures have been suggested besides P(s) 
as providing alternative descriptions of the MIT depend- 
ing on which theoretical and numerical method is being 
used pl| , p8[ . Of particular interest is the so-called number- 
level variance £2, which is a measure of the global spectral 
rigidity |)42"||. It is defined as 



£ 2 (AE) = l[n(AE) - {n(AE))\ 



(4) 



where n denotes the number of levels in a fixed energy in- 
terval AE and () indicates an averaging over disorder. In 
the insulating state £2 = (n>), while it has a logarithmic 
increase £2 oc ln(n) in the metallic state 42 . The behav- 



ior of the number variance at the critical point has been 
conjectured to be Poisson- like |47]], i.e., linear in (n), 



£l 



x(n) 



(5) 



where the level compressibility x is another important pa- 
rameter to characterize the Anderson transition. It is de- 
fined as Esldi 



X 



lim 

(n) — >cx 



lim 

N — >cx 



d£ 2 ({n}) 
d(n) 



(0) 



and takes values < x < 1, being zero in the metallic and 
unity in the insulating state. It is a universal parameter 



and depends only on the spatial dimensionality and on the 
symmetry class [ p0| , |oT| ] . The two limits in Equation (^) do 
not commute. This non-commutativity is attributed to the 



and references 
) of S2 at the 



fractal nature of the critical states (see |4G 
therein). The proposed linear increase (| 
MIT as a function of (n) at the transition has been a 
matter of discussion |58,44|. In general, there is consensus 
that £2 has a quasi- Poisson behavior as in Equation (|^) 
at the MIT 

In this paper we show how £2 and x can be used to- 
gether with FSS to obtain reliable estimates of the critical 
exponent v. We employ various FSS schemes to check the 
accuracy of our results. Our study goes beyond similar pre- 
vious investigations of the level number variance [ j50| due 
to a considerably enhanced accuracy in the scaling data. 
Based on raw spectral data of an anisotropic version of the 
AM, we find that v is consistently larger than 1 in con- 
tradistinction to a recently raised objection |5(| to the FSS 
method of Ref. jL7|, cp. Appendix. The values of v that we 
obtain are in good agreement with the above mentioned 
recent estimates for the isotropic case ]TJ [l9|. |2[ |23|, [37j . 
The mean value of Xc at the MIT is « 0.28 ± 0.06. 



2 The Model Hamiltonian 

We consider the 3D Anderson model of localization de- 
scribed by a Hamiltonian in the lattice site basis as 



H 



^e,-|i)(i| 



K>0'1- 



(7) 



The states \i) are orthonormal and correspond to parti- 
cles located at the A^ 3 sites i — (x,y,z) of a regular cu- 
bic lattice with periodic boundary conditions. The site 
energies Si are taken to be random numbers uniformly 
distributed in the interval [— W/2, W/2]; W defines the 
disorder strength. The hopping integrals ttj are restricted 
to nearest neighbors and depend only on the three spa- 
tial directions. In this paper we consider weakly coupled 
planes defined by t x — t y — 1, t z = 0.1. We emphasize 
that we have chosen the strong anisotropy simply because 
we have the most accurate data (the relative error ranges 
from 0.2 to 0.4%) available for this value from a previous 
study |pl[,|5gj,|59l|. This high accuracy (for spectral data) 
has been achieved by averaging over 10 samples for sys- 
tem size 50 3 and then increasing the number of samples 
up to 699 for system size 13 3 such that always at least 
10 5 eigenenergies have been computed for each N and 
W. Since it was shown in Refs. [£lj|59| numerically that 
the universality class of the model is not changed by the 
anisotropy, we therefore need not generate similarly pre- 
cise data for the isotropic model in order to show scaling 
of £2 and x- 

The Hamiltonian ( |7|) w as diagonalized numerically us- 
ing a Lanczos method |60[. In order to perform any statis- 
tical calculations the eigenspectrum is " unfolded" so that 
the average spacing between adjacent eigenvalues is one. 
Spectra unfolding amounts to a kind of renormalization of 
the eigenvalues in order to extract the universal spectral 
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properties. One way to perform spectral unfolding is to 
subtract the regular part from the integrated density of 
states and consider only the fluctuations 42 1 . This can be 
achieved by different means; however, there is no rigorous 
prescription and the "best" criterion is the insensitivity of 
the final result to the method employed. This criterion is 
fulfilled in the present study. 



3 Finite-Size Scaling 

According to the one-parameter-scaling hypothesis p6[ | , a 
quantity X at different disorders W and energies E scales 
onto a single scaling curve, i.e., 



X(N;W,E) = f{£(W,E)/N), 



(8) 



where the scaling function / is a generalized homogeneous 
function |6^] and £ denotes the correlation length. The 
MIT in the 3D AM is a second-order phase transition and 
as such it is characterized by a divergent correlation length 
with power-law behavior ^ = \W — W c \~ u - The task of 
FSS now is to determine the infinite-size quantities / and 
£(W,E) from finite-size data and to obtain the critical 
exponent v and the critical disorder W c or the critical 
energy E c . 

The essential idea of the FSS procedure of Ref. [0 is 
to construct a family of fit functions which include cor- 
rections to scaling due to an irrelevant scaling variable 
and due to non-linearities of the disorder dependence of 
the scaling variables. The former is only necessary when 
the accuracy of the data allows us to observe systematic 
shifts of the intersection points for different W (or E) and 
N. In all current FSS studies of spectral properties such 
an accuracy has not been reported, only studies using the 
transfer-matrix method allow for an identification of irrel- 
evant variables. 

Following Refs. 



we thus assume a scaling 



form without irrelevant variables to be 



(9) 



where g r is the relevant scaling variable. Taylor expanding 
/ up to order n r , we get 



(10) 



Non-linearities are taken into account by expanding g r in 
terms of u = 1 — W/W c (or u = 1 — E/E c ) up to order m r 



g T (u) = u + h kU k . 



(11) 



k=2 



The fit function is adjusted to the data by choosing the 
orders n T and m r up to which the expansions are carried 
out. 




<n> 



Fig. 1. X , 2((n))/(n) for coupled planes with t z — 0.1 and 
for linear system sizes N = 13,17,21,24,30 and 40 denoted 
by increasing symbol size. The dashed line correspond to the 
GOE. Note that only every 5th symbol is shown. The three 
thick solid lines indicate fits according to Equation (jl^) of data 
for (n) > 15. 



4 Results 

4.1 The spectral rigidity 

In Figure |l| we show the computed £2 data for, e.g., three 
disorders, E G [—4.1,4.1] (50% of the spectrum), and var- 
ious system sizes. In a previous study pl| , we have shown 
that similar level-statistics results can be obtained when 
only 20% of the spectrum close to the band center are 
taken into account p2] . The dependence of W c on energy 
has been considered previously in, e.g. Ref. jfij. A large 
band of states around E = shows the same multifractal 
characteristics as a narrow band |]63| . Thus it is justified 
to take a large part of the spectrum into account when 
computing spectral statistics. It is evident from the figure 
that there is a systematic size dependence as a function of 
disorder. For large disorder W — 12 and upon increasing 
the system size the data approach the insulating (Poisson) 
behavior. Similarly, for small disorder W = 6 the curves 
tend towards the metallic (GOE) behavior. And close to 
the MIT at W c « 8.625, the data for all system sizes col- 
lapse onto a single curve. A similar trend as in Figure |l| 
has been observed for in the four-dimensional isotropic 
AM M. 



4.2 FSS with integrated S 2 data 

In order to perform FSS, we could now use the data of 
Figure [l] and plot them at each value of (n) as a func- 
tion of disorder. However, such an approach is of limited 
usefulness since it is apriori unclear how to weigh data 
from different (n) values. Furthermore, the fluctuations in 
the data lead to rather large error bars in the obtained 
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Fig. 

t» =0.1 for linear system sizes shown. The solid lines are fit 
functions of FSS from Equation (g) with n r = 3 and m r = 1. 



estimates of v and W c . Instead, we define the integrated 
quantity 



V (N,W) = ±- / ° S 2 (L)dL 
L o Jo 



(12) 



with L = (n). This is similar to the FSS analysis of 
/^-statistics |59|. The integral is also considered up to 
Lq = 30 only because the relative error in S 2 becomes 
rather large for larger Lq values and hence the calculation 
is less reliable. It is evident from Figure || that T](N, W) 
shows the desired system-size dependence for various val- 
ues of W exhibiting insulating, critical and metallic behav- 
ior for W larger, close to, and smaller than W c « 8.625. 
In order to obtain £ from finite system-size data we now 
use the FSS procedure of Section 0. For the non-linear fit, 
we used the Levenberg-Marquardt method ]l7| ). In Fig- 
ure U we show that the data from different system sizes 
collapse on two branches corresponding to localized and 
extended behavior. This clearly shows that 77 exhibits one- 
parameter FSS. We then compute the critical exponent v 
and the critical disorder W c for various parameters. The 
results are tabulated in Table [|. The average values are 
v w 1.43 ± 0.13 and W c w 8.62 ± 0.04, respectively. Here 
and in the following, the error intervals are standard er- 
rors, i.e., denoting one standard deviation. 



4.3 FSS with 



We now turn our attention to computing \- First we note 
that S 2 (L)/L (L = (n)) as plotted in Figure |l| is already 
a crude approximation of \- Since there is a systematic 
size dependence, this already indicates that x should obey 
FSS. In order to proceed more accurately, we now fit the 
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Fig. 3. The one-parameter scaling dependence of 77 on £ for 
different system sizes JV and disorders W £ [6, 12]. 

Table 1. Parameters for FSS of 77 and the resulting estimates 
for v and W c . The numbers in the 3rd and 4th column de- 
note orders n r and m r used in the expansions (10) and (11), 
respectively, for which the best fits have been obtained. Q is 
the goodness-of-fit defined as usual by the \ 2 quality of fit Q 
and quoted errors correspond to one standard deviation. 
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^2{L)/L data with an ansatz function containing irrele- 
vant scaling exponents yk, i.e. 



E 2 {N,W,L) 



^a fc (JV, W)L-** 



(13) 



k=l 



up to order m. Thus in the limit L — > 00, the constant 
term will be equal to the desired value of the level com- 
pressibility x- The data used in the fits range from L = 
up to L w 140. Data for larger L values (« 250) was ig- 
nored due to reduced statistical accuracy. In Figure [j], we 
show some typical fits for large system sizes. 

We next perform FSS of \ as explained above. How- 
ever, a non-linear-fit procedure of (increasingly fluctuat- 
ing) large L-data with exponents as fitting parameters is 
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-0.9 1 1 1 1 1 1 1 1 1 

-10 12 3 

Log 10 (N/^) 

Fig. 4. The one-parameter scaling dependence of \, obtained 
from Equation (113) for, e.g., y\ = 1 and y2 — 2, on £ for 
disorders W £ [7jll]> system sizes 24, 30, 40, and 50, with 
n T — 3, m r = 1. The dashed line indicates the value Xc = 
0.27(1) at the MIT obtained from this fit. 



inherently unstable. Thus it is numerically much better 
to fit with fixed exponents. Using such a fit with W — 
7, . . . , 11, N = 24, . . . , 50 for m = 1, we find that \ 
ranges from 0.275 to 0.292 with y\ varying from 0.95 to 
1.3, respectively. Various other combinations for values of 
yi = 0.95, . . . , 1.3 and 2/2 = 1-5,..., 2.5 result in the aver- 
aged FSS estimate Xc — 0.28 ± 0.03. Other values outside 
these ranges do not fit ^(iV, W, L)/L and larger values 
for m do not enhance the quality of the FSS. 

As shown, e.g. for yi — 1 and j/2 = 2 in Figure jij x 
data for different system sizes and W collapse onto a single 
scaling curve with two branches. From this and further 
data for different ranges of W and N, we can roughly 
estimate Xc at the MIT to be 0.28 ± 0.06 as shown in 
Table |[ This value is in good agreement with previously 
obtained estimates (4^,^,p6). We also note that Equation 
(||) with k = l/2xc ~ 1.8 fits the large-s tails of P(s) at 
the MIT reasonably well. 

Unfortunately, the fit in L~ Vk is not good enough to 
reproduce the values for v and W c with the desired high 
accuracy as shown, e.g., in Table ||. This is because the 
data for large L fluctuate much more strongly than at 
small L due to the reduced statistics of such extremely 
large spacings. E.g., even for a system of size 50 3 , we have 
only 250 spacings with L = 250 available when taking 
only the central half of the spectrum into consideration 
to avoid distortions from the localized states in the band 
tails. Furthermore, the usual unfolding procedures rely on 
local spectral interpolations and may no longer work for 
such large spacings. In fact, using an unfolding suitable 
for small L statistics, we could erroneously reduce the es- 
timated value of Xc to w 0.25. 




-0.9 1 1 1 1 1 1 1 1 1 

-10 12 3 



Log 10 (N/^) 

Fig. 5. The one-parameter scaling dependence of x+ > obtained 
from Equation ( |l3| ) for yi — 1 and y% = 2 fitted up to L = 30, 
on £ for disorders W £ [7, 11], system sizes 24, 30, 40, and 50, 
with n r = 3, m T = 1. 

4.4 FSS with x from truncated data. 

In order to suppress the problems with large-L fluctua- 
tions, we have truncated the S^/L data at L = 30 and 
performed the FSS procedure as before with Equation ( jT^ ) 
using fixed y\ — 1 and y% — 2. As shown in Figure ||, the x 
data for different system sizes and different W values col- 
lapse again onto a single scaling curve with two branches. 
Due to the truncation, the estimated values x+ can oniv 
serve as upper limits to the true value % c at the MIT. But 
the resulting values for v and W c are of much better ac- 
curacy and are shown in Table S. For the average critical 
exponent we obtain v w 1.44 ± 0.13 and for the average 
critical disorder W c « 8.66 ± 0.04. 



4.5 FSS with a polynomial fit 

In order to proceed more accurately with the determina- 
tion of v and W c , we now fit E2,{L) data for small L with 
a polynomial in L, i.e., 

Mi^±i aXt+ ± bt{N , w)Lt (14) 

fe=l 

up to order m. We then identify a rough estimate of the 
level compressibility with the linear expansion coefficient 
X p . This implies a systematic shift of x an d the value 
of x P at the MIT will also be shifted towards a larger 
value when compared to Xc- However, v and W c can be 
determined with increased precision: the quality of the fit, 
cp. Figure [| is very good and certainly better than in the 
two previous cases. As a check to the numerical reliability 
of this method, we vary the value of L included in the fit 
function (fj) by fitting the E 2 /L data for various ranges 
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Fig. 6. x„(N, W) obtained by fitting the £2 data with m = 3 
in Equation (|f 4[). System-size dependence is clearly seen. The 
solid lines are fit functions from FSS of Equation (f)|) with n r — 
3 and m r = f . 
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Fig. 7. The one-parameter scaling dependence of Xp on £ for 
different system sizes N and disorders W £ [6, 12]. 



of L = 10, 15, 20, 25 and 30. We find that there is only a 
negligible change in the obtained values of v and W c . 

After FSS, Xp data for all system sizes and all W col- 
lapse onto two curves as shown in Figure Results for 
v and W c for various FSS functions and different m are 
shown in Table |^. For the average critical exponent we ob- 
tain v ps 1.47 ± 0.10 and for the average critical disorder 
W c re 8.56 ±0.05. 

The values of v calculated from the 77- and the three 
X-based approaches are compatible with each other and 
are also comparable to values from other methods P,|67[], 
for instance the transfer-matrix method which gives v ps 
1.62 ± 0.07 f|]]. We can therefore claim that both x and 
77 are good FSS parameters to characterize the MIT. Nev- 



ertheless, a simple fitting procedure in the large L limit, 
although in principle correct, will encounter many numer- 
ical problems. 



5 Conclusion 

States at the 3D MIT are multifractal entities |38], . 
This implies that, while not being extended, their spatial 
structure nevertheless results in a long-r anged, power-law 
overlap of electronic densities in energy |4g,[70|, i.e., 



(|^| 2 |^ fe | 2 ) cx \Ei-E k 



-(l-D 2 /d) 



(15) 



where D 2 is the correlation dimension 1 63 1 and the connec- 
tion Xc = (d— D 2 )/2d has been conjectured [ p0[ . In order 
to describe generic features of such multifractal states, var- 
ious critical random matrix models have been suggested 
and studied |^6],[7l|-f^|, albeit mostly for the unitary class 
of models. Using the above relation of D 2 (see Refs. [|9| 
|63| for numerical estimates at the MIT in anisotropic 3D 
AMs) with Xc we find that the computed value Xc = 
0.28 ± 0.06 is compatible with D 2 « 1.3 ± 0.1, which can 
be calculated easily from the f(a) spectra published in 
Refs. |5^, ^33, 75 1 . On the other hand, it has been shown 
that in the limit of "strong multifractality" , the above 
connection between \c and D 2 no longer holds [^6|. Pre- 
vious estimates of _D(2) in the isotropic 3D AM ran ge fr om 
1.4-1.7 HMfElrfo). Certainly, our multifractal ||, |§ 
is not an infinitely sparse multifractal wave (D(2) = 0)) 
as sometimes expected for the critical ensembles [Em . 

In summary, our results show that x (and nf) can indeed 
be used to compute, with the help of FSS, estimates of 
Xc, W c and v which are in good agreement with transfer- 
matrix and other spectral analysis. We are confident that 
the remaining small difference in values can be further 
shrunk when larger system sizes become available for the 
spectral statistics. 
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A Another FSS procedure 

In a recent communication to the cond-mat archives |p6| , 
the FSS method used in the present paper has been crit- 
icized and the results obtained by various groups |(| 
for the critical exponent v of the localization length at the 
MIT in the 3D AM have been questioned. These claims 
are based on the observation that there still is some dis- 
agreement between analytical, numerical and experimen- 
tal results for the critical exponent H. Ref. p6| proposes 
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Fig. 8. Results for W c and u, for the anisotropic AM |po|,pl[ 
using FSS of Ref. JIt]] for transfer-matrix data (open sym- 
bols) and level-statistics data (filled symbols) for various fit 
functions. The error bars show the 95% confidence intervals. 
The accuracy of the localization length data obtained from 
the transfer matrix method is an order of magnitude higher 
than that of the energy-level-statistics data. The system sizes 
of transfer-matrix data are larger than for level-statistics data, 
giving systematically larger v values for the former. The good- 
ness of a fit is reflected in the size of the symbol. The 2 thick 
error bars mark high quality fits of level-statistics data for large 
system sizes. The gray o and □ and the corresponding dashed 
error bars represent the result of the fit procedure of Ref. 
applied to transfer-matrix data and level-statistics data for the 
anisotropic AM, respectively. The solid line marks the result 
of Ref. 0. 



yet another procedure to deal with corrections to scaling. 
Furthermore, it is hinted that the numerical data support 
v w 1, whereas the present manuscript and recent numer- 
ical papers find v w 1.5 ±0.2 

We have tested the method proposed by Ref. |5(| first 
with transfer-matrix data J2^, £3], |24j with e < 0.1%; we 
find v = 1.75 ± 0.17 for the anisotropic and 1.55 ± 0.04 for 
the random-hopping AM. The FSS of section || gives v = 
1.61±0.07 @ and v = 1.54±0.03 US respectively, for 
the same set of data. Note that the first value (1.75) is so 
high because systematic shifts of W c due to an irrelevant 
scaling variable are not taken into account in | p6| . Using 
for a second test the energy-level-statistics data of the 
present manuscript with e w 1%, we find v = 1.51 ± 
0.25. Last, for artificially generated data with precisely 
known W c — 16.5 and varying v S [0.5, 2.0] the results of 
the method of Ref. |56| are comparable to the results of 
the Kramer/MacKinnon FSS (h]] and slightly less reliable 
than the present FSS as shown in Figure ||. We conclude 
that the method proposed in Ref. J5(| also yields v « 1.58 
and not v » 1 for the MIT of the AM. 
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Table 2. Parameters for FSS of % as in with y± = 1, 
y2 — 2 and the resulting estimates for Xc, v and W c . The 
numbers in the 3rd and 4th column denote orders n r and m r 
used in the expansions (^) and (jTTJ) , respectively, for which 
the best fits have been obtained. The 5th column is for fits 
based on £2 data up to L = 140, the 6 — 8th columns are for 
fits up to L — 30 only. The goodness-of-fit Q |i4j is less than 
0.1 in all cases. 











Lmax = 140 




^max — 30 


w 


N 


n T 


m r 


X 


X+ 


W c 


V 


7--11.0 


13---50 


3 


1 


0.27(9) 


0.31(5) 


8.69(3) 


1.44(7) 


7--11.0 


13- "50 


1 


3 


0.16(3) 


0.31(3) 


8.69(3) 


1.58(8) 


8- -9.25 


13- -50 


3 


1 


0.26(12) 


0.30(5) 


8.69(4) 


1.38(16) 


8- • -9.25 


13- -50 


1 


3 




0.31(5) 


8.68(3) 


1.44(14) 


8- • -9.25 


13- -50 


3 


2 


0.31(1) 


0.33(1) 


8.70(4) 


1.35(15) 


8- • -9.25 


13 -"50 


2 


3 


0.37(2) 


0.37(6) 


8.70(3) 


1.32(16) 


8- "9.25 


24 ■• • 50 


3 


3 


0.27(1) 


0.32(1) 


8.60(3) 


1.11(16) 


7---11.0 


24 ■ • • 50 


3 


1 


0.27(1) 


0.32(1) 


8.60(4) 


1.58(12) 


7---11.0 


24 • ■■ 50 


1 


3 


0.26(2) 


0.32(1) 


8.60(5) 


1.78(10) 


8- • -9.25 


24 ■ ■■ 50 


3 


1 


0.19(13) 


0.32(1) 


8.61(4) 


1.42(17) 


8---9.25 


24 ■• • 50 
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3 


0.36(18) 


0.32(2) 


8.60(4) 


1.45(15) 


7---11.0 


13 "-50 


3 


3 


0.31(1) 


0.33(1) 


8.70(3) 


1.38(10) 


average: 


0.28(6) 


0.32(3) 


8.66(4) 


1.44(13) 
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Table 3. Values of the critical disorder strength W c and the 
critical exponent v values computed from Equation (^) with 
various m values. The numbers in the 4th and 5th column 
denote orders n r and m r used in the expansions (10) and (11), 
respectively for which the best fits have been obtained. The 
goodness-of-fit Q jH^] is larger than 0.99 in all cases. 
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N 


n r 


m r 


w c 
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2 


7--11.0 


13 •• 


50 


3 


1 




L.oO{o ) 
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7---11.0 
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8---9.25 


24 •• 
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3 


1 


8.57(3) 


1.43(6) 
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8- -9.25 


24 •• 
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1.65(7) 
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1.49(5) 
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7--11.0 
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1.42(3) 
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average: 
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24 •• 


50 


3 


1 
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1.56(26) 


4 
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24 •• 
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3 


8.42(11) 


1.68(28) 


average: 


8.48(7) 


1.49(18) 


total average: 


8.56(5) 


1.47(10) 



